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SUMMARY 

In  tnle  paper,  we  wish  to  Introduce  a  new  method  of  treating 
problems  Involving  neutron  multiplication  by  fission,  with  spe¬ 
cial  regard  to  questions  of  critical  mass  and  distribution  of 
neutrons.  Our  results  derive  from  applications  of  the  principle 
of  Invariant  imbedding,  which  we  have  Introduced  and  discussed 
In  two  previous  papers,  the  first  of  which  contains  a  brief 
discussion  of  tne  close  Interrelation  between  the  method  we 
employ  here  and  the  point  of  regeneration  technique  employed  by 


Bellman  and  Harris. 
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ON  THE  PRINCIPLE  OF  INVARIANT  IMBEDDING  AND 
ONE— DIMENSIONAL  NEUTRON  MULTIPLICATION 

By 

Richard  Bellman,  Robert  Kalaba  and  G.  Milton  Wing 
§1 .  Introduction 

In  this  paper,  we  wish  to  Introduce  a  new  method  of  treating 

problems  Involving  neutron  multiplication  by  fission,  with  special 

regard  to  questions  of  critical  mass  and  distribution  of  neutrons. 

Our  results  derive  from  applications  of  the  principle  of  Invariant 

Imbedding,  which  we  have  Introduced  and  discussed  in  two  previous 
1  2 

papers.  *  ’the  first  of  which  contains  a  brief  discussion  of  the 
close  Interrelation  between  the  method  we  employ  here  and  the  point 
of  regeneration  technique  employee  In  Bellman  and  Harris-} 

In  order  to  present  the  essentials  of  the  technique,  as  free 
as  possible  of  analytic  difficulties,  we  shall  consider  first  only 
one— dimensional  flux  of  neutrons.  This  case  Is  of  interest  In  It¬ 
self,  from  both  the  physical  and  mathematical  point  of  view. 

We  shall  consider  first  a  one-velocity  case,  deriving  partial 
differential  equations  for  the  generating  functions  of  the  proba¬ 
bility  distribution  of  "reflected"  end  "transmitted"  neutrons, 

(those  terms  will  be  defined  precisely  below),  and  ordinary  differ¬ 
ential  equations  for  the  moments.  Although  the  moment  equations 
can  be  used  to  determine  critical  length,  an  Inspection  of  the  first 
and  second  moments  shows  that  these  yleld#only  a  very  imperfect 
picture  of  the  actual  distribution  of  neutrons.  Classical  treat¬ 
ment  Is  based  on  expected  values. 

W'e  then  turn  to  a  two— veloclt,  case,  to  illustrate  the 
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appllcablllty  of  our  methods  to  the  more  realistic  situations 
where  fission  and  fission  products  are  energy  dependent,  as 
well  as  situations  where  several  types  of  particles  are  Involved. 

It  will  be  clear  from  our  discussion  and  previous  papers 
how  our  procedure  can  be  modified  to  take  Into  account  absorption 
of  neutrons  and  lnhomogenelty  of  the  medium. 

In  subsequent  papers,  we  shall  first  consider  the  depen¬ 
dence  on  time  of  the  distribution  of  neutrons  and  the  case  cf 
continuous  dependence  on  energy  for  the  one— dimensional  case, 
and  then  turn  to  a  discussion  cf  the  two— d4 mens  Iona  1  plane 
slab,  the  cylindrical  case,  ana  the  spnerlcal  case. 

•  Description  of  the  Process . 

Consider  a  finite  one— dimens  Iona  1  "rod",  tne  Interval 
0,0],  constituted  of  a  material  with  the  property  that  a 
neutron  traversing  an  Infinitesimal  Interval,  [y+dy,y] ,  has,  to 
first  order  terms,  tne  probability  ady  of  splitting  Into  two 
neutrons,  of  the  same  type  as  tne  original,  one  going  to  the 
left  and  one  to  the  right. 

Assuming  a  neutron,  the  "trigger"  neutron,  enters  from 
trie  left,  as  pictured  below,  we  wish  to  determine  the  number  of 
neutrons  emerging  from  the  left,  "reflected"  neutrons,  and  the 
number  of  neutrons  emerging  from  the  right,  " transml tte j"  neutrons 

trigger  neutron  — >  | - 

x  0 


reflected  neutrons  <— 


— »  transmitted  neutrons 


Matnemat lea  1  Formulation  . 
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Let  us  define  the  two  sequences  of  probabilities, 
fpn  (-x  )y  ,  ^  qp(x)£,  as  follows: 

Pn(x)  =■  the  probability  that  n  particles  will  be  reflected 
from  an  interval  of  length  x,  over  all  time,  as  a 
result  of  one  neutron  entering  at  x  at  time  0. 

q^(x)  -  the  probability  that  n  particles  will  be  transmitted 
through  an  Interval  cf  length  x,  over  all  time,  as  a 
result  of  one  neutron  entering  at  x  at  tlm«  0. 

Tne  essence  of  the  technique  of  invariant  imbedding  Is  to 
regard  the  process  described  above  as  a  member  of  a  family  of 
processes  of  tne  same  type.  The  mathematical  statement  of  tnls 
Invariance  will  ylalJ  a  set  of  equations  governing  the  process. 
This  Is  accomplished  readily  In  tnls  case  by  regarding  x  as  a 
parameter  capable  of  assuming  all  positive  values. 

§4.  Basic  Func t lonal  hquat Ions  . 

Let  us  begin  by  deriving  relations  connecting  members  of 
the  sequences  fpn(x+A)y  and.'p^(x)(  for  small  A. 


1 - 1 - - - - - 1 

X  +-A  X  0 

Referring  to  tne  diagram  above,  we  see  that,  to  first  order 


terms,  pQ(x+A)  ( 1-a.A  )pQ  (x  )  . 

To  obtain  the  equations  for  pn(x+A),  we  consider  reflection  of 
n  neutrons  from  ^x^-AjOj  as  composed  of  tne  following  occurrences: 
a.  transmission  of  the  reutron  thru  [x+-A,Xj,  reflection 
of  n  neutrons  from  |_x,Cj,  transmission  of  n  neutrons 
thru  lx  >*■*-&]  • 


Lion 
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b.  transmission  of  the  neutron  through  jx+A,x],  reflection 
of  k  neutrons  from  [x,0],  k«l,2,...,n,  fission  of  one 
of  these  k  neutrons  In  jyc,x+A],  reflection  of  (n— k) 
additional  neutrons  from  [x,0]. 

c.  fission  of  the  neutron  In  [x+A,x] ,  reflection  of  (n— 1) 
neutrons  from  [x,0]. 

This  verbal  statement  Is  equivalent  to  the  mathematical  rela- 

n 


(2) 


Pn(x+A)  -  ( 1— aA )  Lpn(x) (1-aA)  +aA  Z  kpk (x )p  ,  (x  )] 4-aAp  , (x  ) ,  (3) 

k»l 

to  first  order  terms  in  A. 

Letting  A  — >  0,  the  limiting  equations  are 


pQ(x)  -  — ap0 (x ) 


n 


Pn(x)  -  -(n-fl  )apn  (x  )  +apn_x  (x )  +a  Z  (x )pn„k (* ) ,  n  £  1-  (*0 

k  ®  1 


Similarly,  we  obtain  tie  recurrence  relations 


n 


qn  *)  *  -anqn(x)^-a  I  kqk(x)p  k(x),  n»l,2,..., 

k'l 

with  q  x )  -  0 . 


(t) 


§'•  Partial  Differential  Equations  for  Generating  Functions 
Let  ua  define 

v-'x.r)  -  Z  p„(x)rn,  v(x,r)  -  Z  qn(x)rn. 


n 


n-0  11  n»  1 

The  recurrence  relations  of  (4.4)  and  (4. 5)  yield  tne  equations 

ux  -  au(r-l)  ♦  arur(u-l), 

vx  -  ar(u— 1 )vr . 

These  equations  can  be  explicitly  resolved  using  the  method 
of  characteristics.  The  details  are  not  trivial,  however,  we 
shall  discuss  these  matters  In  detail  subsequently. 

§6 .  Critical  Length  and  Moment  Equations . 

The  concept  of  "critical  length"  can  be  introduced  as  the 


(1) 


(2) 


This  Is  equivalent  to  the  fact  that  In  an  Inf lnl teelmal 
Interval  It  Is  sufficient  to  consider  only  first  order  processes. 


smallest  value  of  x  for  which  the  expected  number  of  reflected 
and  transmitted  neutrons  Is  Infinite. 

Setting 

UOO  -  Vr-l’  V<x)  *  (1) 

we  obtain  from  (5*2),  via  differentiation  and  setting  r-1,  the 
equations 

u'(x)  =  a ( 1+U2 ) ,  v'(x)  -  aUV,  (2) 

with  U(o)  =  0,  V(o)  =  1.  Prom  these  equations  we  obtain  U(x)  *  tan  ax, 
V(x)  =■  sec  ax,  showing  that  the  critical  length  le  v/Pq. 

In  a  similar  fashion,  we  find  the  differential  equation  for 
the  second  moment, 

dUp  P 

=  2aU  +  2a  IT  +  3aUU2,  U2 JO)  -  0,  (3) 

which  shews  tnat  U2  =  0  ( 1/ (gob  ax)  ^  )  as  x  ->  r/2a .  It  follows  that 
the  expected  number  of  neutrons  Is  an  unreliable  measure  for  x  near 
T/2a  . 

§ 7 .  Two— Velocity  Case . 

Following  the  usual  approximation  to  the  physical  situation 
where  a  neutron  possesses  a  direction  and  an  energy  varying  contin¬ 
uously  over  certain  limits,  let  us  assume  that  there  are  Just  two 
types  of  neutrons,  "fast"  and  "slow",  and  that  in  the  fission 
process,  either  can  give  rise  to  the  other. 

To  simplify  the  equations,  let  us  assume  that  fast  neutrons 
have  a  probability  apdy  of  splitting  In  {y+dy,yj,  and  that  when 
they  do  spilt  there  Is  a  probability  1/2  of  a  faat  neutron  produced 
going  to  the  right  and  a  slow  neutron  produced  going  to  the  lef t ,  and 
probability  1/2  of  a  slow  neutron  going  to  the  right  and  a  fast 
neutron  to  the  left.  Let  the  same  situation  prevail  with  slow 
neutrons  with  ap  replaced  by  a,.. 
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P ^ mn ( x )  3  th®  probability  that  m  fast  neutrons  and  n  slow  (1) 
neutrons  will  be  reflected  from  an  Interval  of 
length  x,  over  all  time,  as  a  result  of  one  fast 

fc  j 

neutron  entering  at  x  at  time  0,  with  Pv^(x) 
defined  similarly  foY  a  slow  trigger  neutron. 

Then,  as  above,  we  have 

P(mn^+A)  •  h-aA)  [,P(mn<*)(1-aA)m+n+  1  P^  ^x)  [.  ipi^k  ,n-  ^ > 

4  lrJ»  -L*'.  ^  ^  ^  .  (  r  )  I  ' 


4  ka„  -Ip  '  l  # (x  )  1  1 
V  aHrr>-k,n-i'  'jj 

<  aA  -!•  *p(P?  (x)  +  p(S)  ,  (x)] 

8  LKm— l,nx  '  y  m,n-lK  ' J  j' 


with  a  similar  equation  for  p^'(x+A). 

Passing  to  the  limit,  we  obtain  recurrence  relations  which 
permit  us  to  obtain  partial  differential  equations  for  tne  generating 
functions,  u  ^ '  'x  ,  r,  t )  =  Z  p^^jx)rratn  and  u^(x,r,t)  »  Z  p  ^  (x  )  rmtr< . 
Prom  tnene  we  obtain  ordinary  differential  equations  for  the 
moments,  and  thus  can  determine  critical  length.  At  the  present 


we  have  not  been  able  to  solve  these  equations  explicitly,  though 
numerical  solution  would  be  straightforward. 
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